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ABSTRACT 

We explore the possible evolutionary status of the primary component of the bi- 
nary 85 Pegasi, listed as a target for asteroseismic observations by the MOST satellite. 
In spite of the assessed 'subdwarf status, and of the accurate distance determination 
from the Hipparcos data, the uncertainties in the metallicity and age, coupled with the 
uncertainty in the theoretical models, lead to a range of predictions on the oscillation 
frequency spectrum. Nevertheless, the determination of the ratio between the small 
separation in frequency modes, and the large separation as suggested by Roxburgh 
(2004), provides a very good measure of the star age, quite independent of the metal- 
licity in the assumed uncertainty range. In this range, the constraint on the dynamical 
mass and the further constraint provided by the assumption that the maximum age is 
14 Gyr limit the mass of 85 Peg A to the range from 0.75 to 0.82 Mq. This difference 
of a few hundreths of solar masses leads to well detectable differences both in the 
evolutionary stage (age) and in the asteroseismic properties. We show that the age 
determination which will be possible through the asteroseismic measurements for this 
star is independent either from the convection model adopted or from the microscopic 
metal diffusion. The latter conclusion is strengthened by the fact that, although metal 
diffusion is still described in an approximate way, recent observations suggest that the 
real stars suffer a smaller metal sedimentation with respect to the models. 
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1 INTRODUCTION 

The bright visual binary 85 Peg (HD 224930) has been the 
ob ject of numerou s studies since the time of its discovery 
bv lBurhainl lll879l) . Spectroscopy revealed a metal deficient 
object a nd astrometry moderately high space velocity com - 
ponents ('AU ende Prieto et al.l2004HNordstrom et al.i2004 . 
both suggesting a classification as subdwarf and a probable 
old age. Its orbital elements are known with good approxi- 
mation since a long time, and recent determinations, made 
through very precise astrometric. p hotometric and spectro- 
scopic measurements llGriflinll2004l and references therein), 
have not substantially changed what has been known for 
more than one century nor have been able to solve the prob- 
lem of the masses of the stars that constitute 85 Peg. In- 
deed, while it seems evident that, making the reasonable 
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hypothesis of a common origin, the secondary star should 
have a considerably smaller mass than the primary, due to 
the about three magnitudes of difference, the masses deter- 
mined dynamically are similar. This discrep ancy has b een 
ascribed to a third much fainter component |Hal|ll94^. In 
this situation the determination of the two masses has been 
committed up to now to theoretical calibrations by modeling 
stellar parameters (Fernandes et al. 1998, Fernandes, Morel 
& Lebreton 2002) on the basis of observational data and of 
assumptions on age, initial helium content and metallicity. 
This leaves broad intervals of discretionality which, coupled 
with the uncertainty in the theoretical models, lead to rather 
different determinations of mass and consequently of age. 
The determination of the age of subdwarfs sets a lower limit 
to the age of old galactic stars, possibly significant as lower 
limit for the age of the Galaxy. 

Recently, a new tool has become available to infer the 
evolutionary status and the structural properties of a star. 
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asteroseismology. During the last years, numerous attempts 
have been made with the aim of identifying oscillations in 
distant stars and also of modeling the stellar pulsational 
phenomena. It is now clear, indeed, that oscillations have 
several advantages over all the other observables: pulsational 
instability has been detected in stars in all the evolutionary 
stages and of different spectral type; frequencies of oscilla- 
tions can be measured with high accuracy and depend in 
very simple way on the equilibrium structure of the model; 
different modes are spatially confined and probe different 
layers of the stellar interior. 

85 Peg A can be considered as a good candidate for 
asteroseismic studies, since theory predicts a detectable 
surface amplitude of oscillations with a maximum value 
of Voac — 18.8 cm s~^ a ccording to the scaling law of 
iKieldsen fc Beddind (^^S) and a mean value of Vosc — 9 
cm s~^ (Houdek private communication) obtained accord- 
ing to the theoretic al simulation of m odes stochastically 
driven by convection (|lfoudek et al.ll9^fl^ . 85 Peg A, indeed, 
is listed as target ib r asteroseismic observations by MOST 
ijWalker et al.ll200,^ . the first satellite totally dedicated to 
the observation of pulsating stars from space. 

The aim of the present paper is to provide a comprehen- 
sive theoretical study of the system. We will use new the- 
oretical evolution models to find a correspondence between 
the observed data and the derived astrophysical parameters, 
in order to locate 85 Peg A in an evolutionary context. We 
will discuss the possibility of discriminating among its possi- 
ble structure models by means of the oscillation data which 
will be soon available. 

A summary of known data on 85 Peg is presented in 
Section |^ the details of the evolutionary models and their 
result in Sections 1^- 13 the computations of the oscillation 
frequencies and a discussion of the results in Section |S] and 
|7| which is followed by the final conclusions. 



2 A "THORN IN THE SIDE" OF 

EVOLUTIONARY THEORIES: THE VISUAL 
BINARY 85 PEGASI 

The visual and single-lined spectroscopic binary 85 Peg is 
a bright (V=5.75), and nea rby (~ 12 pc ) system with an 
orbital period of 26.31 yr llGrifiinllioM ). Its small angu- 
lar separation (a = (f.'83) and the marked magnitude dif- 
feren ce between the components. Am ~ 3.08 ± 0.29 mag 
in V llten Brummelaar et al.ll200(tl . makes it a difficult tar- 
get both for visual and spectroscopic observations . The pri- 
mary component is a metal deficient G5 subdwarf iFulbrighd 
|2000|), while the most-of-the-time unresolvable secondary 
is supposed to be, on the basis of the Am value, a K6- 
K8 dwaxf. The classification as subdwarf is based both on 
the space velocity component values and on its met allic- 
ity jAllende Prieto et al.l 12004 iNordstrom et al.ll20oj . re- 
spectively [Fe/H ] = -0.88 and -0.80). More recently 
IVan't Veer et al.l ll2005h found a smaller deficit, [M/H] = 
—0.7, which at an y rate places it a t least among old disk 
stars. According to lVan't Veer et al.l (12005), the a-elements 
are enhanced by a factor ~2 with respect to solar scaled 
abundances. 

An excellent summary and a critical analysis of the dy- 
namical data, collected over more than a century, has been 



Table 1. Orbital and astrometric parameters of 85 Peg 



parameter 


value 


error 


ref. 


P (yr) 


26.31 


±0.01 


2 


a 


0'.'83 


itO'.'Ol 


1 


i 


-49° 


±1° 


1 


e 


0.38 


±0.01 


1 


TT (mas) 


82.5 


±0.8 


1 


Ki (km s'^) 


4.49 


±0.05 


2 


7(km s^-"-) 


-36.22 


±0.03 


2 


ai sin i (AU) 


3.68 


±0.04 


2 


f{M) (Mq) 


0.0722 


±0.024 


2 



Reference coding: f P) ISoderhielml il99flh . (2) iGriffinI i2004h 
Table 2. Physical parameters of 85 Peg 



parameter 


primary 


secondary 


ref. 


M/Mq 


0.77 ±0.05 


0.72 ±0.05 


1 


Teff(K) 


5600 ±50 


4200 ±200 


2,3 


Mv 


5.39 ±0.04 


8.47 ±0.29 


4 


L /Lq 


0.617 ±0.02 


0.072 ±0.03 


4 


[M/H] 


-0.7 ±0.1 




2 


log g 


4.6 ±0.1 




2 



Erro rs ar e standard deviations. R efere nce coding: (1 'llOriffinI 
l2004) . f2') IVan't Veer et al.l ^2005^ . (3) iFernandes et al.l i2o'oi) 

(4) present paper 

recently published by iGriffinI ||2004) together with a long 
term radial velocity study, therefore we report here only the 
strictly necessary information and refer to the quoted pa- 
per for a more complete review. The main updated orbital 
parameters from variou s sources are collected in Table [J 

The often quoted llLippincottI Il98ll : IGriffinI 120041) re- 
mark by IStruvd (^5^ on 85 Peg: "this double star is a 
thorn in the side of the present evolutionary theories of 
single stars and of widely separated double stars" essen- 
tially refers to the discrepancy between the reduced mass 
{B — M2/{Mi + A^2)), as derived from astrometry, and 
the value expected on the basis of the magnitude difference. 
The former is typically around 1/2, e.g. t he recent value B = 
0.528 ± 0.034 from Hipparcos results llMartin fc MignardI 
Il998h . and implies a secondary star equally or slightly more 
massive than the primary, the latter is around B = 0.38. 
With a well constrained total mass of 7Vf=1.49±0.1 7Vf0the 
secondary mass varies, then, from 0.77 to 0.55 Mq. 

The value of the reduced mass fr om astrometry is con- 
firmed by the spectroscopic orbit of IGriffinI i2004l) . as his 
value of the mass function f{A4) — 0.0722 A40, when cou- 
pled to the total mass and inclination from visual observa- 
tions, provides B — 0.482 ± 0.019, in good agreement with 
the abovementioned value. A possible and straightforward 
explanation is that the secondary is in its turn a binary sys- 
tem, with an undetected secondary (85 Peg Bb) storing the 
"excess" mass (a few tenths of solar mass) . This hypothesis 
was already put forward bv .Hall (194 8) and reappears in a 
number of subsequent works. 

The atmospheric parameters of 85 Peg have been as 
well the subject of several studies. An accurate determi- 
nation of the primary effective temperature and chemical 
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compositi on, by fittinR of the Balmer lin e wings , was ob- 
tained bv lVan't Veed i200(]t) . as reported in lFernandes et all 
(^200^, hereafter FML02), and has recently been updated 
t|Van't Veer et al. 2005,). According to the latter work, the 
primary effective temperature is Toff= 5600 ±50 K. Other, 
somewhat differ ent, determinations yield also lower values; 
[Fulbrigl^ ll2000t) obtained the lowest one (Toff= 5275 K), 
but the discrepancy is presumably due to the derivation of 
the tem perature from the Fel lin es under the assumption 
of LTE. iThevenin fc Idiard showed that non-LTE ef- 

fects can affect the result in the case of sub-dwarfs and the 
same iFulbriehd l|^00) finds indeed a systematic effect, his 
temperatures being lower than those from other methods. 

The individual absolute magnitudes and relative errors 
were determined by us using the apparent magnitude s and 
relative errors given by Iten Brummelaar et alJ ll2000f) who 
got BVRI differential photometry of the components by 
means of adaptive optics, and derived the individual magni- 
tu des from these and from the composite m agnitudes found 
in lMermilliod. Mermilliod fc Hauckl (Il997l) . We used, more- 
over, the Hipparcos parallax and errors as corrected by 
ISodcrhiclm ( 1999) for binarity effects. The bolomet ric cor- 
rection s were obtained from the best fit coefficients of lFloweil 
l|l996l) . The physical parameters of the binary are collected 
in Table |5| Errors on luminosities do not take into account 
the uncertainty on bolometric correction which should be 
negligible at the primary temperature but could be relevant 
for the secondary component. 

FML02 determined by a least square algorithm the set 
of physical parameters, with their evolutionary stellar code, 
that best fit the available photometric, astrometric and spec- 
troscopic data (under the hypothesis of equal age and chem- 
ical composition of the components). They used, however, as 
input to the procedure, a rather lower value for the reduced 
mass {B = 0.44 vs. B = 0.482 from Griflin 2004). Conse- 
quently they determine a value for the primary mass, 0.88 
Mq, which is significantly larger than Griffin's one. These 
authors derive a best fit age of 9.3 ±0.5 Gyr and a (today's) 
metallicity value [Fe/H] = —0.55. 



3 MODELING 

We use the Aton 3.0 version of the Aton code (Mazz- 
itelli 1979) whose main ph ysical inputs are described in 
I Ventura fc D'Antonal i2005l) . We describe here the proper- 
ties of this version of the code which are relevant to the com- 
putation of the oscillation frequencies: namely the choice of 
opacities, the modeling of convection and diffusion, and the 
computation of an updated equation of state. Only models 
having gray boundary conditions are used in this work. 



3.1 Opacities 

We adopt the [Alexander fc FergusonI (^93) opacities at 
temperatures smaller than 12000 K, and the OPAL opac i- 
ties, in the version documented bv llglesias fc Rogers! il996l) . 
Opacities and thermodynamics are computed in the code ac- 
cording to the prescriptions given in the following Section. 



3.2 Equation of state 

Tables of equation of state (EOS) are computed outside the 
main code. For each metal abundance, represented by the 
mass fraction Z, six tables are built, with hydrogen content 
X from to 1— Z. The input variables are temperature T 
and gas pressure Pgas, and the output quantities are the 
density p, specific heat at constant pressure Cp, the ratio 
7 = Cp/Cv, the adiabatic gradient Vad, Xp s-nd the expo- 
nents Fi, r2 and Fs. Simple functional relations then provide 
any other thermodynamic quantity. Tables are rectangular 
in the plane T-Pgas, in order to simplify interpolation by 
means of third degree bidimensional spline functions. No- 
tice that the regions of very high temperature and very low 
pressure are never used, and are filled up with simple for- 
mulations of the EOS, only to provide more convenient in- 
terpolation. The tables are filled up in three steps: first the 
thermodynamic quantities are computed according to the 
formulation by Stolzmann fc Blocker (1996, 2000) for fully 
ionized gas. This EOS allows to consider explicitly differ- 
ent elements mixtured. ft includes degeneracy (both classic 
and relativistic), coulombian effects, exchange interactions 
and other corrections, and it is the most modern descrip- 
tion available for ionized gas. It holds, however, only for 
high densities. 

The region in which pressure ionization is important 
is overwritten with ISaumon. Chabrier. fc van iHo3(|l99i) 
EOS. In this way we fill the whole region i n which the 
gas is not ionized. However, iSaumon, Chabrier. fc van HornI 
{ 1995) EOS is only given for pure hydrogen and pure helium, 
without metals. Therefore the metallicity Z is "simulated" 
through Y and the different hydrogen - helium mixtures 
are interpolated through the additive volume law. This is 
particularly unsafe where strong configurational effects are 
expected, and the interpolati on error may amount to ~ 20% 
for normal H-rich mixtures jMontalban et alJl2000l) . 

Use of ISaumon. Chabrier^^ar^HornTi^g^ EOS is 
limited to a small region of the T-Pgas plane, by overwriting 
it with t he OPAL 2001^ EOS llRogers. Swenson. fc Iglesiaj 
[1996; Ro gerj200j) . which is given for any Z and ratio H/He. 
This EOS extends also to low pressures. 

After these large tables are written, for a given Z, 
six values of each physical quantities are computed for 
6 different Y. A cubic unidimensional spline provides the 
interpolation for any input value of Y. The six tables 
for H/He and given Z are supplemented by fifteen tables 
He /C/O in which the EOS is d irectly computed according 
to IStolzmann fc Blocker! l!l99d) as the non ionized regions 
are not present in stellar structure following helium igni- 
tion. The interpolation among the fifteen tables is done by 
triangles in the plane C/O, as the stechiometric condition is 
Y=l-Xc-Xo. 

In summary, the structure code computes the EOS 
quantities and the opacities memorizing the spline coeffi- 
cients, and this allows to obtain full continuity both in the 
physical quantities and their first and second derivatives, a 
feature important when we deal with the computation of 
multi-mode pulsations. 



1 the WEB version dated March 2002 is used 
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Figure 1. The location of 85 Peg A (full dot with the small er- 
ror box) and 85 Peg B (box) in the HR diagram, according to 
the data in Table |5] The present solar location is also shown. 
We plot isochrones of t =8 (dash-dotted) and 12 Gyr (dotted) 
lines for metallicities Z=0.004 and Z=0.006. For both sets of 
models Y=0.24. Two evolutionary tracks are also shown (solid 
lines), which bracket the reasonable range of ages for the star: 
0.791 A^0(Z=O.OO4), which has an age of about 8 Gyr at the 
luminosity of 85 Peg A (model 4), and 0.775 A^© (Z=0.006), age 
t ~ 14 Gyr (model 5). 



in the framework of the FST model. This mixing is very 
relevant for hydrogen or helium burning in convective cores, 
but it is not strictly important in this work, as we will be 
deahng only with burning of hydrogen in radiative cores. 

3.5 Numerical resolution 

Particularly to compute the solar structure to be compared 
with the very precise results of helioseismology, we take care 
to have a very good numerical resolution of the structure, 
both in the temporal evolution (about 1000 time- steps are 
taken to reach the solar age) and spatial resolution (about 
2000 mesh points are adopted). For the exploratory models 
relative to 85 Peg A presented in this work, we only employ 
~800 mesh points and broader temporal evolution. In any 
case we took care of having an appropriate distribution of 
meshes both in the central region (to have a careful com- 
putation of nuclear evolution), in the overadiabatic region 
and close to the bottom of the external convective region (to 
have a careful computation of the oscillation frequencies). 

3.6 Use of solar-scaled abundances 

In the code, we did not yet implement the use of a-enhanced 
opacities and EOS. At present, the uncertainties both in 
the metallicity determination and in the metal diffusion (see 
next section) are such that this improvement is not yet nec- 
essary, although it must be planned for the future. For pop- 
ulation II stars, we follow the suggestion by Salaris et al. 
(1993), who have shown that the effect of the increase of a- 
elements abundances in the energy production and opacities 
can be mimicked by increasing the total metal abundance. 



3.3 Diffusion 

The code includes b oth gravitational and the rmal helium 
diffusion according to llben fc McDonaldl (|l985!). Helium dif- 
fusion is a key ingredient in the solar model, necessary to 
adequately reproduce the run of the sound speed as de- 
rived from the ex tremely precise helioseismological data 
jSasu. Pinsonneault fc Bahcall 2000.) . 

3.4 Convection 

Convection in the code can be treated ei- 
ther according to the Full Spectru m Turb ulence 
(FST) model bv ICanuto fc Mazzitellil (ll99lD and 
ICanuto. Goldman fc MazziteUil il^9^ or according to 
the classic mixing length theory (MLT) model (Bohm- 
Vitense 1958). A detailed descriptio n of the convective 
model is given in IVentura et al.l il998l) . The FST model is 
characterized by convective fluxes in which the full eddies 
distribution is accounted for, and by a convective scale 
length defined as the harmonic mean between the distance 
of the layer from the top and the bottom of the convective 
zone. In addition, the distance to the top (and bottom) is 
increased by a fraction PHp, where /3 is calibrated in order 
to obtain the observed solar radius at the solar age (in the 
range /3 = 0.1 — 0.2). Non instantane ous mixing is generally 
used as described bv IVentura et al.l (1998). The mixing is 
treated as a diffusion, with the diffusion velocity obtained 



4 METAL DIFFUSION IN PRESENT 
EVOLUTIONARY CODES 

The current version of the evolutionary code does not in- 
clude metal diffusion, whose effect is however much contro- 
versial. In fact, the same "sign" of metal diffusion is not yet 
clear, particularly in the conte xt of low metallicity stars. The 
problem is widely discussed bv lGratton. Sneden fc Carrettal 
(2004). First of all, the theoretical results differ according 
to the description employed: most results based on the de- 
scription by Thoul, Bahcall, fc Loeb (1994), which is widely 
used, predict for the metallicity of the cluster NGC 6397 
a depletion in [Fe/H] by 0.2 - 0.4 dex (depending on the 
ag e) at the turnoff. On the other hand, the recent models 
bv lRichard et all i2002l) . that take into account the effect 
of partial ionization and radiative accelerations, show that 
whereas some elements like He and Li are expected to be 
depleted, others (like Fe) are expected to be significantly en- 
hanced. From the observational point of view, Gratton et al\ 
i200lfl find that turnoff and subgiant stars in the globular 
cluster NGC 6397 have the same metal abundance within 
0.1 dex. This is neither consistent with model s which predict 
a strong diffusion in main sequence, nor withlRicha rd et all 
I2002J models, as there should be a large difference be- 
tween the turnoff abundance -where the metals have been 
either partially depleted from the whole convective region, or 
enhanced- and the subgiants, where deep convection should 
have restored the initial metallicity. 
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In any case, standard models predict that the effect of 
diffusion at the turnoff of globular cluster stars increases 
when decreasing the metallicity, as the main sequence shifts 
to increasingly hotter temperatures, for which the external 
convective layers are smaller, providing higher diffusion ve- 
locities at their bottom. Therefore, a negligible metal dif- 
fusion is expected for 85 Peg A, in view of the very small 
-if any- effect present in the much more metal poor cluster 
NGC 6397. 

In view of these results, it is difficult to accept such a 
large effect as produced in the models by FML02, which start 
with a quasi-solar metallicity ([Fe/H] = —0.185) to arrive 
at the present metallicity, assumed to be [Fe/H] = —0.55. In 
particular, these models would require that the high velocity 
giants in Hipparcos catalogue have on average metallicities 
larger by a factor two than their turnoff or main sequence 
counterparts, a feature which has not been observed. 

Nevertheless, in order to quantify the effect of the 
metal diffusion, we have also considered models computed by 
employing the Code Liegeois d'Evolution Stellaire (CLES) 
in which the helium and metal diffusion has been imple- 
mented according toErhoul ct al. (1994). Comparisons be- 
tween models calculated with and without metal diffusion 
will be shown in Sec. 7.2. 



5 85 PEG IN THE EVOLUTIONARY 
CONTEXT 

We computed, by using the Aton code, two sets of evo- 
lutionary tracks for 85 Peg A and B with helium mass 
fraction Y = 0.24 and metal mass fraction Z=0.004 and 
Z=0.006 correspondin g, respectively, to solar scaled abun- 
dance [Fe/H] ~ -0.70 IVan't Veer et all2005h and [Fe/H] ~ 
—0.55 (FML02). Considering these two different metallic- 
ities, we also get a first understanding of the uncertainty 
due to the adoption of solar scaled abundances in our code. 
The cross se ction are computed according to the NACRE 
compilation lAngulo et al.lll999ll . All the models apart from 
model 9 in Table |3| adopt the FST treatment of the con- 
vective transport of energy in the interior, with a fine tun- 
ing parameter /3 = 0.17. In order to quantify the effect of 
the treatment of convection on the pulsational properties of 
the star, we also computed a track using the Bohm-Vitense 
mixing length formalism (MLT) with omlt ~ I /Hp = 1.7 
(model 9). The a value has been chosen so that the MLT 
and FST tracks have the same location in the HR diagram. 

According to its space velocity, 85 Peg belongs to a 
group of moderate - high velocity stars, whose age is statis - 
tically limited to be larger than ~8 Gyr iCaloi et al.lll999l) . 
We have then considered 8 t ^ 14 Gyr as possible ages for 
85 Peg. In Tableware collected, for each Z, the mass of the 
model reaching the luminosity observed for 85 Peg A (see 
Table|21l at the assumed ages of about 8, 10, 12 and 14 Gyr. 
Fig. shows the locations of 85 Peg A and B according to 
the data from Table |5| the 8 and 12 Gyr isochrones for the 
different metallicities and two extreme evolutionary tracks 
(model 4 and 5). We notice that all the computed isochrones 
are compatible with the location of 85 Peg B. The 12 Gyr 
isochrones of both metallicities are consistent with the Tofi 
location of 85 Peg A, while the isochrone of Z=0.004 and 
t = 8 Gyr has a too large Tcft. However, we can consider 
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Figure 2. The masses and lu minosities of 85 Peg A and 85 Peg B, 
as derived by iGriffiiJ |200^, are compared with the mass lu- 
minosity relation of the computed isochrones shown in Fig. [Tl 
and with that for t=14 Gyr and Z=0.004 [model 1]. If one ex- 
cludes ages larger than ~ 14 Gyr, the minimum acceptable mass 
is ~ 0.75 Mq. 

as well this last isochrone on the basis of uncertainties on 
temperature. On the observational side, the error of only 
±50 K assigned to the Tejr determination certainly does 
not include consideration of possible systematic errors. From 
the stellar structure modeling side, we do know that precise 
Teff's can not be obtained from first pri nciples. In partic- 
ular, we employ grey atmosphere models: iMontalban et all 
(2001) have shown that non gray models for low metallic- 
ity predict slightly lower Tad than gray models (by ~ 80 K 
at Z=0.0002). Our results adopt the FST convection model, 
but a convective model with a slightly smaller efficiency may 
provide smaller Toff. In summary, the models for Z=0.004 
(the best value for the metallicity determined from spectra) 
and ages smaller than 12 Gyr have also full right to be con- 
sidered in relation to 85 Peg, as we will check by looking at 
the mass luminosity relation. 

The mass luminosity relation provides a better infor- 
mation than the HR diagram, as it is not affected by the 
uncertainty in Teff. The mass luminosity relations of the 
computed isochrones shown in Fig. are given in Fig. |5| 
with also that for t ~ 14 Gyr and Z=0.004. The location 
of 85 Peg A according to the data in Table |2l taking into 
account the uncertainties, is in good agreement with all the 
computed Z=0.004 or Z=0.006 isochrones. Besides, in the 
hypothesis of an age smaller than 14 Gyr, we can put a 
lower limit to the primary mass at 0.75 A4q. In summary, 
we can consider for our analysis all the models at the lumi- 
nosity of 85 Peg A, from 0.75 7\/f0(Z=O.OO4), which has an 
age of ~ 14 Gyr to 0.82 >l0(Z=O.OO6) which has an age of 
~ 8 Gyr. 

Notice the location of 85 Peg B in the figure: indeed its 
dynamical mass is much larger than any acceptable value 
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Table 3. Selected models for 85 Peg A 



mod. 


convection 


M/Mq 


Z 


t(Gyr) 


L/Lq 


R/R0 


Toff 


X, 


1 


FST 


0.750 


0.004 


13.80 


0.616 


0.817 


5661 


0.1670 


2 


FST 


0.763 


0.004 


11.96 


0.617 


0.810 


5689 


0.2103 


3 


FST 


0.778 


0.004 


9.97 


0.619 


0.803 


5716 


0.2669 


4 


FST 


0.791 


0.004 


8.31 


0.618 


0.797 


5737 


0.3240 


5 


FST 


0.775 


0.006 


13.97 


0.619 


0.846 


5572 


0.1653 


6 


FST 


0.788 


0.006 


12.10 


0.618 


0.841 


5587 


0.2100 


7 


FST 


0.803 


0.006 


10.08 


0.618 


0.833 


5615 


0.2675 


8 


FST 


0.820 


0.006 


7.88 


0.616 


0.823 


5644 


0.3443 


9 


MLT 


0.803 


0.006 


10.02 


0.615 


0.833 


5608 


0.2706 



at that Riven luminosity. If we consider the value given by 
iGrifiinI (1200411 , the dynamical mass of the secondary compo- 
nent is ~ 0.2 Mq larger than the possible main sequence 
mass. In the hypothesis that 85 Peg B is in its turn a binary 
system, its primary component, 85 Peg Ba, should have a 
mass of « 0.52 Al© (depending on the age), and its compan- 
ion, 85 Peg Bb, of about 0.2 A^q. Accordingly, the luminos- 
ity of component Bb would be ~ 0.0085 LQand that of the 
primary of 0.064 LQ(as the value in Table |21 would refer to 
the composite luminosity). The binary nature of 85 Peg B 
could be checked by spectroscopic or interferometric obser- 
vations (depending on its unknown period), if an instrument 
capable of separating the components of the visual binary 
is used, such as an adaptive optics telescope. The various 
possibilities in relation to the expected orbital period have 
already been examined in ICriffin ( 2004). 

In the next Sections the models are taken as input to 
derive the oscillation spectrum. 



6 PULSATION ANALYSIS 

We use the adiabatic oscillation code by J0rgen 
Christensen-Dalsgaard available in the WEB to calcu- 
late the p-mode eigenfrequencies. Some analytical and 
numerical aspects of oscillation codes are described in 
IChristensen-Dalsgaard and BerthomierJ il99l|) . 

From the asteroseismic point of view 85 Peg A can 
be classified as a solar-type star, a main-sequence star in 
which oscillations are excited stochastically by vigorous 
near-surface convection in a broad spectrum of low ampli- 
tude p-modes, as in the Sun. The interpretation of the pul- 
sational properties of a solar-type star can be o btained by 
considering the asymptotic theory jTassouJll98'ol) . This pre- 
dicts that the oscillation frequencies Vni of acoustic modes, 
characterized by radial order n at harmonic degree Z, should 
satisfy the following approximation: 



1 1\ 

n+-+a+- + e,ii 

2 4/ 



(1) 



where a is a function of the frequency determined by the 
properties of the surface layers and tni is a small correction 
which depends on the conditions in the stellar core. A is the 
inverse of the sound travel time across the stellar diameter: 

-or")"' 

where c is the local speed of sound at radius r and R is 



the photospheric stellar radius. The property expressed by 
Eq. (0 may provide almost immediate asteroseismic infer- 
ences on stellar parameters and constraints on theoretical 
models for a variety of solar-like stars in a wide range of 
evolutionary stages. For a given I, the acoustic spectra show 
a series of equally spaced peaks between p modes of same 
degree and adjacent n, whose frequency separation repre- 
sents the so called large separation which is approximately 
equivalent to A: 



(3) 



The spectra are characterized by another series of peaks, 
whose narrow separation is dii+2, known as small separa- 
tion: 



d, 



I 1+2 



= Vnl - Vn-1 1 + 2 = (4/ -|- 6)Do 



where 



Do 



47r2!/„ 



=GR) 
R 



dc dr 
dr r 



(4) 



(5) 



A, and hence the general spectrum of acoustic modes, scales 
approximately as t he squa re root of the mean density, that 
is, A oc \/p = yj M/ R^ . On the other hand, the small 
frequency separations are sensitive to the chemical compo- 
sition gradient in central regions of the star and, hence, 
to its evolutionary state. Thus, the determination of both 
large and small frequency separation. A; and dii+2, pro- 
vides measures of the mass and of the age of the star (e. g. 
IChristensen-Dalsgaardl[l988ll . 



6.1 Sub-surface effects and the ratio of small to 
large separation 

The seismic analysis of observed acoustic frequencies is an 
extremely powerful tool for the investigation of the inter- 
nal structure of the stars, but the use of large and small 
separations can be misleading if not accurately considered. 
Theoretical pulsation frequencies, essential for asteroseismic 
investigation, are calculated on theoretical models of stars 
which are inevitably affected by errors. In particular, the 
structure of the near-surface regions of stars, which strongly 
affects the large separation, is quite uncertain. In fact, there 
are still substantial ambiguities in the theoretical descrip- 
tion of: i) the convective flux, ii) the mechanisms of exci- 
tation and damping of the oscillations, iii) the equation of 
state to describe the thermodynamic properties of the stellar 
structure, and iv) the non-adiabatic effects. Limiting the 
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investigation to the use of small separation does not solve 
the problem, since the small separation, which is determined 
principally by the conditions in the core, also retains some 
sensitivity to the mean density and to the detailed properties 
of the stellar envelope. 

To solve these pro blems, the use of a new se i smic i ndica- 
tor was introduced bv [Roxburgh fc Vorontsovl i2003ll . who 
compared the pulsating properties of several models with 
exactly the same interior structure, but with different outer 
envelopes. They showed that the ratios of small to large sep- 
arations, 



n = di 1+2/ Ai 



(6) 



are independent of the structure of the outer layers and, 
hence, can be used as diagnostic of the interior of stars. 

6.2 Echelle diagram 

In order to study the large and small separations of the ob- 
served and computed oscillation frequency spectra, it is com- 
mon to use the echelle diagram iGiec. Fossat fc Pomerant j 
^S^). In this diagram each frequency is expressed in terms 
of an integer multiple of A, according to 



- kA + Vnl 



(7) 



where uq is an arbitrary reference frequency, k is an integer 
and Vni is a residual frequency which lies in the range — 
A. Plotting uq -\- kA against Uni one expects, according to 
the asymptotic theory (Eq. , that computed frequencies 
fall in parallel columns, one for each harmonic degree. The 
distance between consecutive frequencies of same harmonic 
degree, belonging to the same column, represents the large 
separation; modes of even degree (0 and 2) and modes of odd 
degree (1 and 3) fall in two pairs of closely spaced columns, 
the distance between the columns in each pair representing 
the small separation. 



7 PREDICTION OF THE OSCILLATION 
FREQUENCIES OF 85 PEG A 

We computed eigenfrequencies for all the selected structure 
models described in Table |3 and consistent with the ob- 
served basic parameters of the primary component of 85 Peg. 
Since the oscillation modes which can be detected in stars 
are generally characterized by low harmonic degree (owing 
to the point-like character of the sources), we limited our 
analysis to the calculation of modes with harmonic degrees 
I = 0,1,2,3. Theoretical calculations on the selected mod- 
els predict that 85 Peg A is pervaded by acoustic modes 
with frequencies in the range 1000 — 5000 /^Ifz for harmonic 
degrees Z = — 3. The oscillation spectrum is character- 
ized by a large separation for the radial modes Aq of about 
150- 170mHz (Fig.EJ. 

As predicted by theory, the value of the large separation 
decreases as the mean density decreases. For fixed metallic- 
ity, the large separation increases with the mass. In fact, 
the large separation strongly depends on the condition on 
the surface, and hence on the value of the metallicity. This 
means that the value of the mass cannot be estimated from 
the observed large separation alone. 

It is very instructive to consider the behaviour of the 
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Figure 3. Open symbols refer to models with Z=0.004, filled 
symbols to Z=0.006, equal symbols correspond to models of simi- 
lar age. Upper figure: large separation Aq for I = as function of 
the frequency calculated for selected models described in Table ISl 
Lower panel: ratio of small (io2 to large separations Aq between 
I = and I = 2 for eight models of Table 1!^ 



large separation as a function of frequency (Fig. |^. Con- 
sistently with the solar case, the theoretical A; shows an 
evident oscillatory dependence on frequency. Such quasi- 
periodic behaviour of the frequencies rises from discontinu- 
ity's regions localized at certain acoustic depth in the struc- 
ture. This peculiar property of the oscillation frequencies 
can be used to infer the helium abundance in the stellar 
envelope or to determine the properties a nd the location of 
the base of the con v ective envelope (e.elBa su et alJIgOOj: 
iMazumdar fc Antial 1200 ll : iMonteiro fc ThompsonTligod : 
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Figure 4. The top figure shows the hydrogen mass fraction X 
along the structure of selected models for increasing age. The 
bottom figure shows that the two models having age close to 
10 Gyr have an almost equal hydrogen profile in spite of the 
different metallicity. 



Figure 5. Large separation Aq for 1 = (top figure) and ratio of 
small to large separation (bottom) for model 7 and 9 of Table 3, 
having same mass and radius and nearly the same age of 10 Gyr. 
The models differ in the modelling of the external convective lay- 
ers: EST for model 7 and MLT for model 9. 



Monteir o Christensen-Dalsgaard fc ThompsonI Il994 Il998l 
2Q0Q , 20(m 

Fig. El employs the param eter ro = do2/Ao devised by 
iRoxburgh fc Vorontsovl ll2003l) . in order to subtract the con- 
tribution of the surface layers, yielding a diagnostic of the 
stellar interior alone. In fact, the difference in cio2/Ao among 
the models is easily understood in terms of the different 
internal structure: independently of the metallicity of the 
models, it is the core hydrogen content -and thus the age 
reached at the luminosity of 85 Peg A- which determines 
do2/Ao. Thus, models in the same evolutive state have the 
same value of ri . The plot of hydrogen mass fraction X inside 
some of the examined structures is shown in Fig. 2] 

It is clear that while measurements of small separation 
(or ratio r;) will be useful to define the evolutive state of the 
star, the large separation will be employed to identify, among 
structures of same age, elements abundances and mass of the 
model which better reproduces the observations. 



T.l MLT versus FST models 

While we adopt the FST model as a standard for the overa- 
diabatic convection, we have made a comparison between 
model 7, computed with FST, and model 9, computed with 
the MLT. We adopt a = I / Hp=l.l , which provides a track 
having the same Tcff location of the FST track. The com- 
parison is made in Fig.|^ there is a difference of about l/xHz 
in the large separations for Z = 0, while the ratios ro of the 
two models are very similar. Thus we conclude that the con- 
vection treatment does not affect the age determination by 
the seismic analysis. 



7.2 The metal diffusion 

The ATON code does not include metal diffusion, so we 
adopt an entirely different code to have an idea of the in- 
fluence of this physical input on the structure. The stellar 
structure and evolution code CLES adopts OPAL 2001 EOS; 
OPAL96 opacity tables plus Alexander & Ferguson (1994) at 
T < 6000 K; nuclear reaction rates from Caughlan & Fowler 
(1988); MLT convection treatment; atmospheric boundary 
conditions given by Kurucz (1998) at T = Teff. Microscopic 
diffusion o f all the elements i s computed by using the sub- 
routine bv lThoul et"ai] Jl994) . A detailed comparison of the 
results of the CLES and ATON codes is out of the purpose 
of this paper, although it is in our research projects. We 
list in Table 4 the physical parameters of two models in 
the error box of 85 Peg A. We choose a large age (12 Gyr) 
to make the comparison, in order to maximize the effect of 
diffusion. The initial metallicity is Z=0.006, corresponding 
to [Z/X]o=-0.555. Model 11 includes only helium diffusion'^, 
while model 10 includes both helium and metal diffusion. At 
the age of 12 Gyr [Z/X]s, the vale at the surface of model 10 
has been reduced by ~33%^, a relatively mild effect which 

^ Notice that the CLES model has Xc=0.20 in the center, while 
the corresponding ATON model 6 has Xc=0.21, in spite of being 
slightly older. The treatment of helium diff usion differs in the 
two codes. The description in the ATON code Jiben &: McDonaldl 
[h)853 provides smaller diffusion velocities and this leads to this 
difference in Xc after 12 Gyr. 

^ Notice that this decrease in [Z/X]s in the surface layers is due 
both to the decrease in Z and to the increase of X due to the 
helium diffusion (this can be seen, e.g. in Fig. 5 and from the 
change of [Z/XJj, also in model 11 which has constant Z). 
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Figure 6. Comparison between the CLES models showing the 
effect of metal diffusion on the large separation (top) and on the 
ratio of small to large separation (bottom). 



confirms the prediction made in Section 4 that metal dif- 
fusion on a star having the Teft of 85 Peg A should not 
be dramatic. We show in Fig. |S] the oscillation properties 
of two models. We can see that the large separations differ 
by ~ 2/iHz, but the ratio ro is again very similar. Thus we 
conclude that the uncertainty in the computation of metal 
diffusion in the case of 85 Peg A is totally contained into 
the uncertainty in metallicity, and does not affect the age 
determination. 



7.3 Echelle diagrams 

In Fig. Q we compare the echelle diagram for model 4 of 
85 Peg A with that of a standard solar model having Y=0.28 
and Z=0.02. The theoretical echelle diagram of 85 Peg A 
has been obtained by using the value of A = 170 fiHz, and 
uq — 870 /iHz. The lower panel of Fig. Q obtained for the 
Sun, is drawn for both the theoretical and the observed 
frequencies by using the value of vq = 690 /^Hz as a ref- 
erence frequency and a large separation A — 136 /^Ifz cor- 
responding to the observed average which can be calculated 
on MDI/SOHO frequencies (Schou 1998). Like in the Sun, 
the theoretical echelle diagram of 85 Peg A shows a regular 
behaviour of the frequencies, indicating no departures from 
the simple asymptotic expression. In addition, it is likely 
that the observed oscillation spectrum of 85 Peg A will not 
show presence of g modes nor of mixed modes. In fact, our 
theoretical calculation shows that g modes are well confined 
in the core and with an amplitude which is too shallow to 
be detected at the surface. It is important to notice that the 
observed solar data at high frequencies deviates from the 
model predictions. This is a consequence of a not fully ade- 
quate treatment of the surface layers (see Sec. 16.111 . There- 
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Figure 7. Echelle diagrams respectively for Model 4 (upper 
panel) of 85 Peg A and a model of the Sun together with the obser- 
vational results (lower panel). The filled symbols show observed 
frequencies, while the open symbols represent the computed fre- 
quencies. Circles are used for modes with I = 0, triangles for / = 1, 
squares for 1 = 2, and diamonds for I = 3. Theoretical frequencies 
of 85 Peg A are plotted with A = 170 /xHz, and uq = 870 /^Hz. 
Theoretical and observed frequencies of the Sun are plotted with 
A = 136 AiHz and i/q = 690 ^lUz. 



fore this behaviour is also expected for 85 Peg A which is, 
like the Sun, a low effective temperature star. 



8 CONCLUSIONS 

In this paper we have addressed the problem of studying the 
structural properties and identifying the evolutionary state 
of 85 Peg A. The location in the H-R diagram indicates that 
this star is in the main-sequence phase of evolution, with a 
mass in the range A1=0.75-0.82 Aleand with an age which 
can vary from 8 to 14 Gyr. These values have been obtained 
by assuming the most recent observational parameters of 
temperature and luminosity, but the classical method of fit- 
ting the stellar modeling parameters to the observational 
data does not allow to reduce the uncertainty on age and 
mass. 
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Table 4. Selected models for 85 Peg A obtained with CLES code with and without metal diffusion 



mod. 


convection 


metal diff. 


M/Mq 


[ Z/X]o 


[Z/X]s 


t(Gyr) 


L/Lq 


R/Rq 


Teff 




10 


MLT 


YES 


0.788 


-0.555 


-0.678 


11.90 


0.615 


0.843 


5572 


0.2035 


11 


MLT 


NO 


0.788 


-0.555 


-0.583 


11.96 


0.611 


0.851 


5538 


0.2036 



Here we have shown that the problem of determining 
the evolutionary state of 85 Peg A will be solved once ob- 
servations of the oscillation spectrum will be available for 
this star. Theoretical calculation of models and p mode os- 
cillations have shown that 85 Peg A is a solar-like star, 
which is expected to oscillate in the range of frequencies 
1000 — 5000 ^Hz. Our results show that its oscillation spec- 
trum seem to be characterized by a large separation for I — 
of about 150—170 /iHz and a small separation between I — 
and I = 2 which can vary from about 2 to 6 fiRz, according 
to the evolutionary state of the star. 

Indeed a difference of about 2 Gyr corresponds, in the 
asymptotic regime and hence at high frequencies, to a differ- 
ence in Ao of about 3 fiHz and in the ratio do2 / Ao of about 
0.007 (see Fig.OJ. It can be concluded that a determination 
of frequencies with an uncertainty ^ 0.8 ^Hz, which corre- 
sponds to an average error of ±0.005 in the ratio do2/Ao, 
is sufficient to determine the age with an uncertainty better 
than 2 Gyr. 

Our work has examined the role of the most important 
physical inputs (metallicity, convection and diffusion pro- 
cesses) in the framework of the inputs of our evolutionary 
code. We have studied the influence of metal diffusion with 
the help of the CLES evolutionary code, which differs from 
ours in some physical inputs (e.g. nuclear reaction rates, 
treatment of helium diffusion, atmospheric boundary condi- 
tions, convection modelling). We therefore did not perform 
an accurate comparison of the codes, but simply have shown 
that the inclusion of metal diffusion has not an important 
impact on the age deter mination. Anyway, the change in 
metallicity implied by the lThoul et alJ (Il994ll description of 
metal diffusion is relatively modest (~ 33%) in 12 Gyr, while 
FML02 find a reduction of metallicity by more than a factor 
2 in ~ 9 Gyr. This difference is taken by us as an indication 
that really much has to be done in the understanding of this 
issue. 

We believe that this star should be considered as pri- 
mary target for space or ground-based observations. It is 
clear indeed that only accurate observations will allow in- 
vestigation of the properties of the interior of this star and 
discrimination of its accurate evolutionary stage. The oscil- 
lation spectrum of 85 Peg A will reveal how old is the star 
in terms of core hydrogen consumption, thanks to the very 
good mass constraint coming from the astrometric and spec- 
troscopic observations. In this respect, it appears feasible to 
put a lower limit to the age of the Galaxy based on the 
asteroseismology of 85 Peg A. 
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